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In most studies of quantum channels, it is assumed that the errors in each use of the channel 
are independent. However, recent investigations of the effect of memory or correlations in error 
have led to speculation that non-analytic behaviour may occur in the capacity. Motivated by these 
observations, we connect the study of channels with correlated error to the study of many-body 
systems. This enables us to use many-body theory to solve some interesting models of correlated 
error. These models can display non-analyticities analogous to quantum phase transitions. 



An important problem in quantum information is the 
determination of the channel capacity of noisy quantum 
channels [ij. In a typical scenario, we wish to send in- 
formation over many uses of a noisy quantum channel. 
At a cost of lowering the information content per parti- 
cle, it can be shown that quantum error correction can 
essentially eliminate all errors. This leads to quantum 
analogues of the channel capacity - the optimal rate at 
which information may be transferred with vanishing er- 
ror in the limit of many uses. In this work we will be 
concerned with the capacity for transmission of quantum 
information which we denote Q{£) for a channel £. 

Most works on quantum channels assume that the 
noise is independent in successive transmissions. How- 
ever this is never exactly true and in many realistic sys- 
tems there can be correlations in the noise. Such channels 
are termed 'memory channels' and their capacity may 
be significantly affected by memory effects (see [Sl for 
a recent experiment). Consequently, quantum memory 
channels have received considerable attention recently 
(see e.g. 0, [1, H, 0] and references therein). In most 
models that have been considered explicitly (e.g. [H]), 
the correlations in the noise are modelled by a small 
number of memoryvarameters. Initial investigations of 
prototype models have suggested that the capac- 
ity may undergo sharp, non-analytic, changes at certain 
values of these parameters. These investigations have 
not been conclusive, however, as with current techniques 
these models cannot be analysed in the relevant setting 
of a very large (in fact infinite) number of channel uses. 

The main aim of this work is to show that for a vari- 
ety of interesting models this obstacle can be overcome 
by relating the study of memory channels to the study 
of many-body physics. An advantage of the framework 
that wc introduce is that the criticality of the underlying 
many-body systems can manifest itself as non-analytic 
behaviour in the channel capacity of the corresponding 
memory channels, thus proving the conjectured existence 
of such effects in memory channels. This suggests, as 
demonstrated here, that the methods of many-body the- 
ory may be used in the study of memory channels to ob- 
tain new results. In the first part of the paper we discuss 
a general framework for the channels that we consider, 



discussing how a capacity formula may be derived. We 
then provide an explicit bound applicable to random uni- 
tary [20[ channels in terms of a thermodynamic quantity. 
In the case of dephasing channels that are random appli- 
cations of orthogonal dephasing unitarics this is exact. 

Mapping many-body systems to memory channels. - 
A standard way of describing noise is to assume that 
each transmitted 'system' particle interacts via a uni- 
tary U with its own environment. In order to introduce 
memory effects we will modify this approach by assert- 
ing that the environment particles are initially prepared 
in the thermal or ground state of an interacting many- 
body Hamiltonian. The spatial correlations in the envi- 
ronment then lead to correlations in the noise. We will 
assume that once the environment has been defined by 
the parent Hamiltonian, no further dynamics occur other 
than the system-environment interaction. It is important 
to note that many of the noise models considered in the 
literature 0, H, S, 0] can be reexpressed in precisely this 
way. 

For the case of a memoryless quantum channel it has 
been shown that the quantum capacity given by 

Q{£) = hm -^^-^ (1) 

n^oo 77, 

where £ is the channel acting on a single transmission, 
I[£) :=supp [S[£{p)) - S{I ® £{\^){'il}\))] is the coherent 
information of the quantum channel £ where S denotes 
the von-Neumann entropy, p is a state, and |V')('0l is a pu- 
rification of p. Finally, f*^" represents the uncorrelated 
channel that acts on n inputs. For channels with corre- 
lations, the channel on n inputs £n differs from £f^, and 
one has to describe the memory channel by a sequence 
of channels {£n}, describing the action of the channel 
for each number of inputs n. Given that eq. ^ is the 
memoryless quantum capacity, one may anticipate that 

Q({f„}) := lim (2) 

n — ^oo 77 

is the quantum capacity of a memory channel. While this 
will certainly not be true in general, eq ([2|) is always an 
upper bound on the channel capacity and one can derive 
conditions for equality which are often satisfied: 
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Condition for eg. (0j - This section may be omitted 
by readers not concerned with detailed proofs. We out- 
Hne the derivation of two conditions on the many-body 
system which are sufficient to demonstrate that eq. ([5]) is 
the quantum capacity. These conditions are satisfied by a 
variety of many-body systems, including matrix product 
states , and quasi- free bosonic systems 12 1 . The con- 
ditions are derived as follows. We consider a translation 
invariant chain of length N, split into v = N / {I + s) sec- 
tions, each consisting of one live block of length / and one 
spacer block of length s := 61 <^ I . We define two related 
channels: the Live channel and the Product ch.a.nne\. The 



Live channel is Si, 



= A^ tI^nv{UipL^L2....L^ ® 

where A is Alice's input, U is the interaction between 
chain and input, the ii, . . . , L„ label the live blocks, and 
the environment is traced out. The Product channel is 
defined as Sproduct := A ^ tr,nv{U{{p^j^)^^ ® A)W}, 
where p^v reduced state of an individual live 

block. The next three steps labelled (A),(B) and (C) 
construct an analogue of the arguments made in [7| 
for forgetful channels. (A) Showing that product clas- 
sical codes are good classical codes for the Live Chan- 
nel. Given an achievable rate R for the Product chan- 
nel, then Ve > : 3N^ such that for n > chan- 
nel uses there is a n/-qubit code {pi\i^..M and decoding 
measurement {Mi}i....,y for v = [2"'^J such that Vi : 
tT:{£product{pi)Mi\ > 1 — £. Tlic samc procedure used for 
the Live channel then yields, via the triangular inequal- 
ity, tr{fz„,(p,OM,} > 1 - e - \\\pl^l,....l^ - {p'nT^'Wi- 
This leads to our first condition on the many-body sys- 
tem: it turns out that one can ensure that this error 
vanishes, for a choice v = l^ [2lj, provided that it can be 
shown that \\pL^L,....L^ - (p'iv)®1li < Cvl^cxp{-Fs) 
for positive constants C,E,F (see [12|). Hence if this 



condition holds the classical codes for product channels 
are also good codes for the Live correlated channel. (B) 
Computing achievable classical rates. The product chan- 
nel with live block length / and a total number of spins 
N = v{l + s) ~ 1^{1 + 6) has a Holevo quantity given 
by x(^Ar) = x(trenv{C^((p^) «) •)C/^}), where the • acts 
as place holder for the channel input and where f j^- de- 
notes the effect of the full channel upon a contiguous 
subset of j < X of the input spins. As in Q we must 
now understand when this expression converges to the 
regularized Holevo bound of the full memory channel as 
I oo. Suppose that we have a spin chain of total 
length / -I- A{1) where A(/) > is any function such that 
limj_,oo ~ 0. Using subadditivity and the Araki- 

Lieb inequality we find x(.£j+A) — xi^i-i-A) — 2Alog(d). 
Now wc need to show under which conditions this re- 
mains true if the subset of / spins is drawn from a much 
longer chain of length N = 1^{1 -\- S). For any input u; 
to the live block in question the output states will dif- 
fer by at most ||trcnv{t^[w ® {p\+a - P\<i(i+s))^^^}\\i - 
\\U[L,®ip\_,^-p\s^,_^,^)]W\\i < P(/,A) := ||pU-p5vIIi- 



Combining this with x{£i+a) — xi^i-i-A) — 2Alog((i) 
and the Fannes inequality yields limj^oo x('^i6(i+5))/' ^ 
Xoo — hm;^tx) 2P(/, A(Z)) log(d). This gives the second 
condition on the many-body system: if we can pick A(Z) 
such that limj^oo A(0// and limi^oo P{1, ^il)) = 0, 
then the regularized Holevo quantity is the correct clas- 
sical capacity. (C) Coherentification. The final step is to 
argue that the above arguments for classical coding can 
be 'cohcrentificd' Q into a quantum code. This analysis 
does not give new conditions on the many-body system 
and can be conducted as in (see [H for details). 

Explicit computation of capacities - Even if the many- 
body system can be well understood, the explicit compu- 
tation of the capacity may still be difficult, as it depends 
upon the interaction of each system with its environment 
U. Judicious choice of U will allow us to obtain analyt- 
ically solvable models. To this end we will choose U to 
be of the form of a controUcd-PHASE gate, denoted by 
Uz, where the environmental particles act as controls. 
In this case it becomes possible to write down explicit 
formulae for eq. ^ in terms of properties of the many- 
body environment that share a close relationship with 
thermodynamical quantities. For all other random uni- 
tary noise the approach leads to lower bounds, although 
they are not always exact - see the concluding section 
for discussion. For d-dimensional systems the controUed- 
PHASE gate is defined as = Y1=i \k) {k\(Si Z {k) where 
Z{k) :— exp(27rifcr/d)|r)(r|, and the first tensor fac- 
tor acts on the environment. This interaction leads to 
channels that are probabilistic applications of Zik) uni- 
taries on the system particles, with the (correlated) prob- 
abilities determined by the diagonal elements of the en- 
vironment state. Now eq. ([2]) can be written as 



Q({£„}) =logd- 



lim 

n—*oo 



S{T)ia,g{penv)) 



(3) 



where Diag(pe„t,) is the state obtained by eliminating all 
off-diagonal elements of the state of the environment in 
the computational basis. Hence computing the capacity 
of our channel {£n} reduces to computing the regularised 
diagonal entropy of the environment. Although this is 
unlikely to be generally computable, it is amenable to a 
great deal of analysis using many-body theory. 

Proof sketch of equation (Q) The proof utilises the 
Choi-Jamiolkowski representation of the quantum chan- 
nel. Given any quantum operation £ acting upon a d- 
level quantum system, one may form the correspond- 
ing Choi-Jamiolkowski (CJ) state J{£) = /®5(|+)(+|), 
where |+) = X]i=i dN)- '^be proof of equation ^ 
follows from three steps: (1) we argue that for the kinds 
of channels we consider, a copy of J{£) allows one to 
physically implement the channel exactly, (2) for any 
channel that can be implemented using J{£), we argue 
that the quantum channel capacity Q{£) of the channel 
equals Di{J{£)), the one-way distillable entanglement of 
the state J{£), (3) wc then use known results on Di. 
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Step (1) - For simplicity we describe the argument for 
channels £ that are random applications of Pauli rota- 
tions on a single qubit. The argument generalises easily 
to Pauli channels on many qudits. Suppose that you have 
J{E) and you want to implement one action of E upon 
an input state p. This can be achieved by teleporting p 
through J{£). This will leave you with a state £{aipaj), 
with the Pauli error at depending upon the outcome of 
the teleportation measurement. As f is a random Pauli 
channel we can now "undo" the error by applying the 
inverse of ai. Hence we have ai£{i7ipa\)i7i = £{p) and 
we have implemented one action of £. 

Step (2) ~ Our aim is to show that for channels that 
may be physically implemented using J{£), the 1-way 
distillable entanglement of J{£), Di{J{£)), is equiva- 
lent to Q{£). In [13 it was essentially shown that 
Q{£) > Di[J{£)). For the converse inequahty consider 
the specific protocol: (a) Alice prepares many perfect 
EPR pairs and encodes one half according to the code 
that achieves the quantum capacity Q{£). (b) She tele- 
ports the encoded qubits through the copies of J{£). in- 
forming Bob of the outcome so that he can undo the effect 
of the Paulis. (c) This effectively implements the channel 
£ between Alice and Bob. (d) Bob decodes the optimal 
code, thereby sharing perfect EPR pairs with Alice, at 
the rate determined hy Q{£). As this is a specific one-way 
distillation protocol, this means that Q{£) < Di{J{£)). 
These arguments extend straightforwardly to any chan- 
nel that is a mixture of Paulis on many qudits. 

Step (3) - The CJ states of our channel are so-called 
maximally correlated state, for which the distillable en- 
tanglement is known and is given by the Hashing bound 
[III Di{J{£)) = S{tvB{J{£)}) - S{J{£)) where S is the 
von-Neumann entropy. Hence for such channels £ this 
expression is also the single copy coherent information. 
In our cases we are interested in the regularised value 
of this quantity, which is given by equation ([3]). This 
expression for the coherent information has an interest- 
ing interpretation for the dephasing interactions that we 
consider - it represents the classical information lost to 
the environment that is needed to correct the errors [l3| . 

The simplicity of eq. Q enables one to immediately 
write down many models for which eq. ([2]) can both be 
calculated, and also represents the quantum capacity. 

Classical Environments - We discuss briefly two cases. 
If the environment consists of classical systems described 
by a classical Markov Chain then in a large number of 
cases eq. ([3]) can be written explicitly with a simple ex- 
pression that represents the entropy rate of the Markov 
chain 1J| . Related results on Markov chain models have 



been obtained using different methods in [15| . In general 
a classical environment is represented by a diagonal state 
and the second term of eq. ([3]) is precisely the entropy. 



Hence in this case the capacity becomes 

Q{{£n}) = 1 - log2(e) (l - /3 J.) lim - lnZ„ (4) 

\ op J n^oo n 

where Z„ is the partition function for n environment 
spins, 13 = l/ksT, and the log2(e) converts from nats 
to bits. Equation (j4]) shows that we can now use results 
from classical statistical physics to compute the capacity 
- any classical spin-chain models with sufficiently decay- 
ing correlations that can be solved exactly will lead to 
memory channels that can be 'solved exactly'. 

Quantum Environments - For quantum environments 
eq. ([3]) represents the entropy that results when every 
environment qubit is completely dephased. Although this 
quantity is not standard in statistical physics, we expect 
that it may be amenable to the techniques of many-body 
theory. Here we provide support for this claim by solving 
analytically a class of quantum environments inspired by 
recent work on Matrix Product States (MPS) [l6|. For 
MPS states the two conditions required to prove equation 
(2) are satisfied except at transition points fl^ . 

In [17| it was shown that there are Hamiltonians 



that exhibit quantum phase transitions and have ground 
states that are MPS involving only rank-1 matrices. We 
will now show that for MPS involving rank-1 matrices 
a full analytical treatment of the quantum channel ca- 
pacity of the associated memory channel becomes pos- 
sible. To this end we demonstrate that the diagonal 
elements of such rank-1 MPS are given by the proba- 
bilities of microstates in related classical Ising chains. 
For simplicity we will focus on a translationally invari- 
ant MPS for a ID system of 2-level particles with peri- 
odic boundary conditions. Generalisation to other rank- 
1 MPS is straightforward. Such an environment state is 
characterized by two matrices Qq and Qi and is given by 

IV') = Z]ii...,ijv • ■ ■QiN}Vi ■ --iN)- The unnormal- 

ized state resulting from dephasing each qubit is 



JV 



■ iN){ii 



.in\ 



(5) 



k=l 



Relabeling the matrices Ai = Qi Q*, the diagonal 
elements in the computational basis are of the form 
trjj^j, Afc}. As the Ai are both rank-1 with unique non- 
zero eigenvalues a^, the normalised matrices Ai = Ai/ai 
satisfy Af = Ai. Using this idempotency it is easy to 
show that if \ii . . .i^) has I occurrences of and N — I 
occurrences of 1, and K boundaries between blocks of Os 
and blocks of Is, then the corresponding diagonal element 
of p will be p{l,n~l,K) = {a^b^-%^{{AoAi)^} /C{N), 
where C{N) is a normalisation factor. Noting that 
AqAi is also rank-1, denote its non-zero eigenvalue by 
c. Hence the diagonal elements are p{l,n — l,K) = 
a^b^~^c^ /C{N). Hence for channels described by rank-1 
MPS states, a,6, c are the only relevant parameters and 
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we may choose any matrices with those parameters. We 
choose the matrices of a classical (i.e. diagonal) chain: 



Ao 




(6) 



Here we have assumed that c > (which is guaranteed for 
iV > 5 as otherwise the state can become non-positive). 
It is then easy to check that these matrices have the cor- 
rect values of a, b, c, as required. We choose this form 
because the matrices are essentially the top row and bot- 
tom row of a transfer matrix 0] corresponding to classical 
Ising chain. Roughly speaking, the parameter c encodes 
the coupling between adjacent antiparallel spins, and the 
a,b encode the couphng between adjacent parallel spins. 

This connection implies that the limit in eq. 131) can 
be computed easily using well known methods [8[. Fig . 
[T] shows the result for a Hamiltonian presented in [17| 
for which the ground state is known to be a rank-1 MPS 
possessing a non-standard 'phase transition' at g = 
[l8|, at which some correlation functions are continuous 
but non-differentiable, while the ground state energy is 
actually analytic Fig- [1] shows that this is mirrored 



Diverging gradient 




g=0 



FIG. 1: A sketch of the capacity for an environment that is 
the ground state of the Hamiltonian 2{g^ — l)cr*(T*^^ — 



{l+gYal + ig-ir 



The plot's symmetry is 



expected as the channel is invariant under g —g. However, 
near the 'phase transition' p = 0, the gradient diverges. 

in the non-analyticity of the channel capacity. 

Generalisations and Future work - It is important to 
know whether our approach could prove useful for other 
interactions. Some generalisations are immediate. For 
instance, given any channels that are probabilistic appli- 
cations of unitaries, expression ^ can easily be shown to 
be an explicit lower bound to the coherent information, 
and hence if the environment state has sufhciently decay- 
ing correlations it will also be a lower bound to channel 



capacity. It is likely that any channel whose capacity can 
be bounded by such simple entropic expression will bene- 
fit from similar insights. In the long term one might spec- 
ulate that there may be a deeper explanation for these 
connections - not in terms of entropic expressions appear- 
ing in both fields, but in terms of a link between coding 
and many-body physics. 
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